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Abstract 

The forward-backward charge asymmetry for the process qq —)• Zj'^* —)• is 

measured as a function of the invariant mass of the dimuon system. Measurements 
are performed using proton proton collision data collected with the LHCb detec¬ 
tor at y/s = 7 and 8TeV, corresponding to integrated luminosities of lfb“^ and 
2fb“^ respectively. Within the Standard Model the results constrain the effective 
electroweak mixing angle to be 

sin^^^ = 0.23142 ± 0.00073 ± 0.00052 ± 0.00056, 

where the first uncertainty is statistical, the second systematic and the third theo¬ 
retical. This result is in agreement with the current world average, and is one of the 
most precise determinations at hadron colliders to date. 
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1 Introduction 


In the Standard Model (SM), the Z boson conplings differ for left- and right-handed 
fermions. The difference leads to an asymmetry in the angnlar distribution of positively 
and negatively charged leptons produced in Z boson decays. This asymmetry depends 
on the weak mixing angle ( 6 *w) between the neutral states associated to the U(l) and 
SU(2) gauge groups, i.e. the relative coupling strengths between the photon and the Z 
boson. In order to compare directly with previous experimental determinations, a scheme 
is adopted in which the higher order corrections to the Z boson couplings are absorbed 
in effective couplings. The resulting effective parameter sin^^^ is defined as a function 
of the ratio of the vector and the axial-vector effective couplings of the Z boson to the 
fermions involved [^, and is proportional to sin^ 6 'w 

Dehning Q* as the polar angle of the negatively charged lepton in the Collins-Soper 
frame, in which the direction of the 2 ;-axis is aligned with the difference of the incoming 
proton momentum vectors in the dimuon rest frame, the differential cross section in the 
SM at leading order is 


dcr 


d cos Q* 


= kl(l + cos"r) + 5cosr. 


Here A and B are coefficients that depend on the dimuon invariant mass, mainly because of 
interference between Z and 7 * contributions, the colour charge of the quarks and the vector 
and axial-vector couplings. The parameter 5 is a function of sin^^w a.nd is proportional 
to the forward-backward asymmetry Hfb, which is given by 


Hfb = 


Nf-Ne 
Afp + -/Ve 


where A^f represents the number of forward decays (cos0* > 0 ) and A^b the number of 
backward decays (cos 6 ** < 0). The Collins-Soper frame is used because it minimises the 
impact of the transverse momentum of the incoming quarks on the identihcation of forward 
and backward decays. 

In this paper the asymmetry of the angular distribution of muons in Z —)■ decay^ 

is measured using proton proton collision data collected by the LHCb experiment at 
centre-of-mass energies of ^/s = 7 and 8 TeV, corresponding to an integrated luminosity 
of lfb“^ and 2fb“^ respectively. The asymmetry as a function of the dimuon invariant 
mass is used to determine sin^^^. 

Comparisons of the determinations of the weak mixing angle from processes with 
different initial and final state fermions provide a test of the universality of the fermion 
to Z couplings. The most accurate measurement of sin^^^ at the LEP experiments 
was obtained from the forward-backward asymmetry in b quark final states [^, and at 
the SLD experiment by measuring the left-right asymmetry with polarised electrons [^. 
Determinations of sin^^^ have also been obtained in hadronic production processes with 

^In the following Z is used to denote the Z/ 7 * contributions. 
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leptonic final states at the CDF and DO experiments at the Tevatron and ATLAS 
and CMS experiments at the LHC 

Measurements of Afb can be related to sin^6*^ when the momentum direction of the 
initial quark and antiquark are known. At the LHC the momentum direction of the initial- 
state quark is not known, diluting the ability to determine sin^^^ from Afb- However, 
since at LHC the dominant production process is uu, dd —?• Z, the main contribution 
originates from a collision of a valence quark with high momentum and a sea antiquark 
with lower momentum, and so the Z boson tends to be boosted along the direction of 
the quark. This is particularly true in the forward region where the Z boson has large 
longitudinal momentum. Consequently, the sensitivity of Afb to sin^0^ is greater at 
large rapidities of the Z boson. Using simulated samples, it is found that in the LHCb 
acceptance the assignment of forward and backward decays is correct in 90% of the time. 

The layout of this paper is as follows. Section describes the LHCb detector and the 
data samples used in the analysis. The candidate selection and background determination 
are described in Sec. 0 In Sec. I^the Afb measurements are presented and in Sec. [^the 
measurements are compared to next-to-leading order (NLO) theoretical predictions within 
the same kinematic region, and a value of sin^6*^ is determined. 


2 Detector and datasets 

The LHCb detector is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < t] < 5. The detector includes a high-precision tracking system consisting of a 
silicon-strip vertex detector surrounding the pp interaction region, a large-area silicon-strip 
detector located upstream of a dipole magnet with a bending power of about 4 Tm, and 
three stations of silicon-strip detectors and straw drift tubes placed downstream of the 
magnet. The magnet polarity can be reversed, so that detector-induced asymmetries can 
be studied and corrected for in the analyses. The tracking system provides a measurement 
of momentum, p, of charged particles with a relative uncertainty that varies from 0.5% 
at low momentum to 1.0% at 200GeVj^ The minimum distance of a track to a primary 
vertex, the impact parameter, is measured with a resolution of (15 + 29 /pt) h-ni, where px 
is the component of the momentum transverse to the beam, in GeV. Different types of 
charged hadrons are distinguished using information from two ring-imaging Cherenkov 
detectors. Photons, electrons and hadrons are identified by a calorimeter system consisting 
of scintillating-pad and preshower detectors, an electromagnetic calorimeter and a hadronic 
calorimeter. Muons are identihed by a system composed of alternating layers of iron and 
multiwire proportional chambers. The online event selection is performed by a trigger, that 
consists of a hardware stage, based on information from the calorimeter and muon systems, 
followed by a software stage, which applies a full event reconstruction. For this analysis, 
candidates are triggered by at least one muon with momentum larger than 10 GeV. 

Simulated samples are used to estimate the shapes of the invariant mass distributions 
for the simulated signal sample and some of the background sources. The signal sample is 

^Units where the speed of light is set to unity are used throughout this paper. 
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also used to correct the data for reconstruction and detector effects. In the simulation, pp 
collisions are generated using Pythia 8 10,11 with a specihc LHCb configuration 12 


Decays of hadronic particles are described by EvtGen 13 , in which hnal-state radiation 
(FSR) is generated using Photos I^. The interaction of the generated particles with the 


detector, and its response, are implemented using the Geant4 toolkit (^, as described 
in Ref. 16 . 

To simulate Z —)■ /i’*'/!” decays with different values of sin^0^, the next-to-leading order 
generator Powheg-Box fl7|, interfaced to Pythia for the parton showering, is used. 


Additional simulated samples are generated without parton showering using Powheg-Box, 
Herwig 18 and Fewz 19 and are used to evaluate theoretical uncertainties. Predictions 


are also obtained using Fewz at NLO and are used to provide an alternative calculation 
of ApB to compare to those computed by the Powheg-Box generator. The same parton 
density function (PDF) is used for both generators. 


3 Event selection 

Dimuon candidates, consisting of two oppositely charged muons, are selected using the 
same criteria as in Ref. (^, but with an extended mass range. The two muons must 
be within 2.0 < rj < 4.5, have good quality track hts, a transverse momentum greater 
than 20 GeV and must combine to an invariant mass within 60 < < 160 GeV. These 

requirements dehne the kinematic region of this measurement. 

The purity of the candidate sample is determined by estimating the contribution from 
background sources using a combination of simulation and data-driven techniques, and is 
found to be greater than 99%. The total yield, reconstructed dimuon invariant mass and 
Afb are determined for each source of background. 

The largest background contributions come from semileptonic heavy-flavour decays and 
events containing misidentihed hadrons, where hadrons punching through the calorimeters 
to the muon stations are identihed as muons, or hadrons have decayed in flight. Both 
contributions are estimated using data-driven techniques. Two heavy-flavour enriched 
samples are selected by widening the mass window and requiring evidence that (i) the 
muons are produced away from the primary vertex, or (ii) that the muons are surrounded 
by hadronic activity. These two samples are combined to estimate both the shape of 
the reconstructed distribution and the total number of events for the heavy-flavour 
background source. The misidentihed hadron contribution is estimated by using a sample 
of same-sign muon events. The Z —)■ ^ tt, single top and W~^W~ —)■ 

background sources are estimated using simulation. The total background contribution is 
largest at low invariant mass. The charge asymmetry of each background component is 
consistent with zero over the whole mass range. The distribution of the dimuon invariant 
mass is shown for data and all background sources in Fig. 
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Figure 1: Dimuon invariant mass distribution for data, simulated signal and background processes 
for (left) a/s = 7 TeV and (right) ^/s = 8 TeV. 


4 Forward-backward asymmetry measurements 


The forward-backward asymmetry is calculated from the selected dimuon candidates. 
Corrections are applied to account for efficiencies, biases in the reconstructed momenta of 
the muons and differences in resolution between simulation and data. Previous studies 20 - 


22 have observed a dependence of trigger, track reconstruction and muon identffication 


efficiencies on muon pseudorapidity. To account for this dependence, correction factors 


are evaluated from data using a tag-and-probe method 20 and applied to the measured 
forward-backward asymmetry. 

The momentum measurement of a muon is sensitive both to uncertainties in the 
detector alignment and the magnetic held scale. The magnetic held scale has been 
calibrated using dimuon and other resonances at low mass and is known to a precision of 
0.04% [^. Low-mass resonances have also been used as input to the detector alignment, 
leading to a well-understood momentum calibration for low-momentum tracks [9,23,24 


However, studies for the analysis presented here, have revealed a small, but appreciable, 
dependence of the position of the Z resonance peak on muon kinematics. This ehect 
can be attributed to residual detector misalignment. The corresponding muon curvature 
bias can be ehectively parameterised in bins of the azimuthal angle of the muon about 
the beam axis. The parameters are determined using the diherence between the Z mass 
peak in data and simulation. The procedure is applied separately to data collected at 
fs = 1 and 8 TeV, and for the two magnet polarities. The results are consistent with 


those presented in Ref. 21 in which a slightly different method was used. 


To compare with theory predictions, the data are unfolded for acceptance and resolution 


effects. A Bayesian unfolding technique 25 is applied to the reconstructed dimuon invariant 


mass distribution [^. The unfolding algorithm is trained on simulation by comparing 
the generated invariant mass to that after reconstruction. The simulation is corrected 
to have the same resolution as observed in data. Finally, the data are corrected for 
background by subtracting the distribution for each background source determined as 
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Table 1: Weighted average of the absolute systematic uncertainties for ^fBj for different sources, 
given separately for y/s = 7 and 8 TeV. 


Source of uncertainty 

75 = 7TeV 

7s = 8 TeV 

curvature/momentum scale 

0.0102 

0.0050 

data/simulation mass resolution 

0.0032 

0.0025 

unfolding parameter 

0.0033 

0.0009 

unfolding bias 

0.0025 

0.0025 


described in Sec. No correction is applied to the measured values of Afb to account for 
the dilution due to imperfect knowledge of the initial quark direction, or to remove FSR 
effects. Instead, they are compared to predictions made within the same kinematic region 
and including FSR, as described in Sec. 

The following systematic uncertainties are considered when determining Rfb- The 
systematic uncertainty associated with the curvature correction is evaluated by varying 
these parameters within their uncertainty. The uncertainties on the calibration factors 
are dependent on the sample size, and are therefore larger for the ^/s = 7 TeV dataset. 
This is the largest source of systematic uncertainty. An uncertainty of ±0.04% is used 
for the momentum scale, determined from measurements of the magnetic held |^. The 
bias in the unfolding procedure is determined from simulation by comparing unfolded 
samples with the generated true distribution. An additional uncertainty to account 
for the dependence on the number of iterations used in the training of the unfolding 
algorithm is determined. This variation has a larger effect in regions where fewer events 
are simulated. The asymmetry of each background source does not vary signihcantly over 
the invariant mass range. An uncertainty of 10% is assigned to the background asymmetry, 
that covers the huctuations observed in Afb for each background source. The effect of 
the uncertainties in the efficiency corrections applied to the data is found to be negligible. 
The systematic uncertainties are determined separately for each bin of invariant mass and 
for both datasets. Their average values are summarized in Table [Tj 

The resulting measurements of Afb for a/s = 7 and 8 TeV data as a function of 
are shown in Fig. |^and tabulated in Tables and 

5 Determination of sin^^|y^ 

The forward-backward asymmetry as a function of the dimuon invariant mass is compared 
with several sets of SM predictions generated with different values of sin^0^, denoted as 
The predictions are generated using Powheg-Box with sin^0^ values ranging from 
0.22 to 0.24 for ^/s = 7 and 8 TeV, and the Z boson mass {Mz) and the electromagnetic 
coupling constant (oem) fixed to the world average values [^. The PDF set from 


5 





1 0.6 
0.5 1- 
0.4 I- 
0.3 1- 
0.2 1 - 
0.1 l- 
0 l- 
- 0.1 1 - 
- 0.2 
- 0 . 


1-1-1-1-1-1-1-1-1-1-1-1-1-1-1-1-1-1-1“ 


LHCb 

p^>20 GeV 
2.0<t|<4.5 


J < 




—(— data <ls = 7 TeV 
, POWHEG + PYTHIA 
(sin^ejf = 0.2315) 


%0 80 100 120 140 160 

mpp [GeV] 


0.6 

0.5 

0.4 

0.3 

0.2 

0.1 

0 

- 0.1 

- 0.2 

- 0 . 


1-1-1-1-1-1-1-1-1-1-1-1-1-1-1-1-1-1-1“ 


T 

LHCb 

p^>20 GeV 
2.0<r|<4.5 


data ^ = 8 TeV 
POWHEG + PYTHIA 
(sin^e^' = 0.2315) 


%0 80 100 120 140 160 

mpp [GeV] 


Figure 2: The measurements of as a function of the dimuon invariant mass for data compared 
to SM predictions for (left) ^/s = 7 TeV and (right) ^/s = 8 TeV. The SM predictions are calculated 
using PowHEG interfaced with Pythia for parton showering with the world average value for 
~ 0.2315 27 . The data include both statistical and systematic uncertainties, and the 


sm 


SM predictions include the theoretical uncertainties described in Sec. 


NNPDF g 


.JB with the strong coupling constant as{Mz) = 0.118, was used when generating 
the samples. 

Theoretical uncertainties associated with the distributions are taken into account 
when determining sin^0^. They arise from the underlying PDF, the choice of renormali¬ 
sation and factorisation scales, the value of as used, and the FSR calculation. Each of 
these uncertainties, referred to collectively as theoretical uncertainties, are obtained from 
simulation. The same uncertainty is assigned to App'^ at both a/s = 7 and 8 TeV. 

To estimate the theoretical uncertainty of the PDF set, one hundred replica samples 
are produced, each with a unique PDF set provided by NNPDF 29 . The value of 


RpB^ is calculated as a function of for each of these replicas, and the corresponding 


confidence level interval determined. The size of this uncertainty is larger than the 
difference observed using CTIO as an alternative PDF parameterisation. 

Uncertainty in the PDFs affects Rps'^ in a way that is correlated across all dimuon 
invariant mass bins. The same systematic uncertainty is applied for both collision energies 
and is therefore fully correlated for the two samples. 

The uncertainty due to the choice of renormalisation and factorisation scales is studied 

The uncertainty in the sin^0^ determination 


31 


by varying them by a factor of 0.5 and 2 
due to the uncertainty in ag is estimated by studying the impact of a variation of ± 0.002 
when generating samples using Powheg-Box. This covers the current uncertainty on 
For both the a., and scale uncertainties the hnal uncertainty is estimated by htting 


a 


27 


a constant across the mass range to the maximum and minimum deviations in Rpp^ to 
minimise the effect of statistical fluctuations in the samples. 

The uncertainty due to the implementation of FSR is treated as a theoretical un- 


3NNPDF 2.3 QCD -h QED NLO. 
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Table 2: Values for ApB with the statistical and positive and negative systematic uncertainties 
for y/s = 7 TeV data. The theoretical uncertainties presented in this table, corresponding to the 
PDF, scale and FSR uncertainties described in Sec. affect only the predictions of ApB and the 
sin^0^ determination, and do not apply to the uncertainties on the measured ^fb- 



(GeV) 

Afb 

stat. 

syst. -f 

syst. — 

theoretical 

60 

- 72 

-0.248 

0.018 

0.011 

0.006 

0.025 

72 

-81 

-0.144 

0.015 

0.006 

0.003 

0.011 

81 

-86 

-0.078 

0.013 

0.005 

0.005 

0.005 

86 

-88 

-0.017 

0.013 

0.007 

0.009 

0.005 

88 

-89 

0.016 

0.013 

0.012 

0.008 

0.005 

89 

-90 

0.023 

0.010 

0.006 

0.006 

0.005 

90 

-91 

0.033 

0.008 

0.007 

0.004 

0.005 

91 

-92 

0.047 

0.008 

0.009 

0.002 

0.005 

92 

-93 

0.082 

0.010 

0.004 

0.010 

0.006 

93 

- 94 

0.127 

0.014 

0.004 

0.016 

0.009 

94 

-98 

0.175 

0.012 

0.003 

0.014 

0.009 

98 

- 120 

0.259 

0.015 

0.007 

0.006 

0.014 

120 

- 160 

0.451 

0.037 

0.004 

0.017 

0.026 


certainty. It is obtained by comparing the from three different generators, Fewz, 
Herwig++ and Powheg-Box+Pythia, before and after FSR. To be consistent with 
the Powheg-Box sample, the Fewz generator is confignred at NLO and electroweak 
corrections are not inclnded. The maximum and minimum difference is then determined 
and used to estimate the systematic uncertainty associated with FSR. The average size 
of the separate theoretical uncertainties is summarised in Table and the combined 
uncertainties as a function of invariant mass are given in Tables and 

The shapes from Powheg-Box were cross-checked using the Fewz generator 
at the same value of sin^^^ and the differences were found to be negligible. 

The agreement between data and prediction is quantihed by a value dehned as the 
square of the difference between the measured RpB and Rpp'^ divided by the quadratic 
sum of the statistical, systematic and theoretical uncertainties, taking into account the 
correlations in the uncertainties between the mass bins. A quadratic fu n ction is fitted 
to the values of each set of fits as a function of sin^0^. The result is shown in Fig. 
The value of sin^0^ at which takes its minimum is quoted as the final result for the 
sin^0^ determination. The interval in sin^ 6 *^ corresponding to a variation of one unit 
in is quoted as the uncertainty. The observed minimum values for the x^/ndf of the 
fit are 0.59 and 0.58, for the 7 and 8 TeV samples, respectively. The minima correspond 
to sin^^^ = 0.23219 ± 0.00148 and sin^ 6 *^ = 0.23074 ± 0.00123 respectively. Results are 
cross-checked using a set of pseudoexperiments with the same statistics and background 
fractions as those in data. The values of App from the pseudoexperiments are fitted to the 
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Table 3: Values for ApB with the statistical and positive and negative systematic uncertainties 
for y/s = 8 TeV data. The theoretical uncertainties presented in this table, corresponding to the 
PDF, scale and FSR uncertainties described in Sec. affect only the predictions of ApB and the 
sin^0^ determination, and do not apply to the uncertainties on the measured ^fb- 



(GeV) 

ApB 

stat. 

syst. + 

syst. — 

theoretical 

60 

- 72 

-0.217 

0.014 

0.015 

0.014 

0.025 

72 

-81 

-0.154 

0.012 

0.004 

0.004 

0.011 

81 

-86 

-0.046 

0.010 

0.003 

0.002 

0.005 

86 

-88 

-0.004 

0.010 

0.003 

0.004 

0.005 

88 

-89 

- 0.002 

0.011 

0.003 

0.007 

0.005 

89 

-90 

0.016 

0.008 

0.006 

0.002 

0.005 

90 

-91 

0.040 

0.006 

0.005 

0.003 

0.005 

91 

-92 

0.053 

0.006 

0.004 

0.002 

0.005 

92 

-93 

0.075 

0.008 

0.004 

0.006 

0.006 

93 

- 94 

0.104 

0.011 

0.003 

0.006 

0.009 

94 

-98 

0.166 

0.010 

0.005 

0.006 

0.009 

98 

- 120 

0.280 

0.012 

0.006 

0.002 

0.014 

120 

- 160 

0.412 

0.027 

0.005 

0.009 

0.026 


Table 4: Weighted average of the absolute systematic uncertainties for ^pB^’ different 

sources of theoretical uncertainty. The value quoted for the PDF uncertainty corresponds to the 
68% confidence range, while for the others the maximum and minimum shifts are given. The 
correlations among the invariant mass bins are not taken into account. 


Uncertainty 

average A 

PDF 

0.0062 

scale 

0.0040 

Otg 

0.0030 

FSR 

0.0016 


prediction, and the spread of the measnred sin^^^ valnes agrees with the nncertainties 
in the valnes of the 7 and 8 TeV samples. A combination of these results, taking into 
account the correlation between systematic uncertainties for each centre-of-mass energy as 
well as the invariant mass bins, is obtained by calculating the full covariance matrix for 
the statistical, systematic and theoretical uncertainties. This yields 

sin^^^ = 0.23142 ± 0.00073 ± 0.00052 ± 0.00056, 
where the hrst uncertainty is statistical, the second systematic and the third theoretical. 
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0.228 0.229 0.23 0.231 0.232 0.233 0.234 0.235 

sin^e^ 

Figure 3: Difference between the ^ and the minimum obtained by comparing the final 
ApBi'iTT'iJ.fj.) measurements in data to calculated using values of sin^ 0 ^ ranging from 0.22 

to 0.24, indicated by the crosses on the plot. A quadratic fit is used to determine the minimum 
value for sin^ 0 ^ and the corresponding uncertainty, and is shown for the different centre-of-mass 
energies and the combination. The black dashed horizontal line corresponds to one unit of 
from the minimum and the intersecting sin^ 0 ^ for the combination are indicated by the vertical 
red dashed lines. 


A comparison between the sin^^^ result obtained here and those from other experiments 
is shown in Fig. The LHCb result agrees well with the world average and is one of the 
most precise measurements from hadron colliders. 


6 Conclusions 

The forward-backward asymmetry for the process gg —)• Z —)• jjA as a function of the 
dimuon invariant mass is measured with the LHCb detector using proton proton collision 
data collected at centre-of-mass energies of y/s = 7 and 8 TeV. The measurements are 
performed in the Collins-Soper frame, using muons with pp > 20 GeV and 2.0 < p < 4.5 
with a combined invariant mass 60 < < 160 GeV. The forward-backward asymmetry 

for each invariant mass bin is measured, together with the statistical and experimental 
uncertainties. The measurements at each centre-of-mass energy are used to determine a 
value for sin^^^, by comparing to SM predictions that include FSR. The best fit values 
obtained are sin^^^ = 0.23219 ± 0.00148 and sin^ 6 *^ = 0.23074 ± 0.00123 for the two 
samples at ^/s = 7 and 8 TeV respectively. This leads to the combined result 

sin20^ = 0.23142 ± 0.00073 ± 0.00052 ± 0.00056, 
where the first uncertainty is statistical, the second systematic and the third theoretical. 
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Figure 4: A comparison of the sin^0^ measurement at LHCb and other experiments. The 
combined LEP and SLD measurement is indicated by the vertical yellow band. 


The measurement of sin^6*^ presented here agrees with previous measurements. The 
uncertainty from the PDF is the dominant theoretical uncertainty. Further high precision 
measurements at the LHC are expected to provide additional constraints in the forward 
region and reduce this uncertainty. As the size of the data sample increases, it will become 
possible to perform a measurement of ApB double-differentially in dimuon invariant mass 
and rapidity. Such an approach will allow the analysis to take further advantage of the 
increased sensitivity of Apu to sin^^^ in the most forward region. 
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